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We investigate a celebrated problem of the one dimensional tight binding model in the presence 
of disorder leading to Anderson localization from a novel perspective. A binary disorder is assumed 
to be created by immobile, heavy particles that affect the motion of the lighter, mobile species 
in the limit of no interaction between mobile particles. Fast, periodic modulations of interspecies 
interactions allow us to produce an effective model with small diagonal and large off-diagonal disorder 
previously unexplored in cold atom experiments. We present an expression for an approximate 
Anderson localization length and verify the existence of the well known, extended resonant mode. 

We also analyze the influence of nonzero next-nearest neighbor hopping terms. We point out that 
periodic modulation of interaction allows disorder to work as a tunable band-pass filter for momenta. 


I. INTRODUCTION 

Anderson localization (AL) in disordered systems has 
fascinated and stimulated physicists for more than 50 
years [1, 2]. For one-dimensional (ID) systems, as dis¬ 
cussed below, a particle propagating with momentum k, 
in a disordered medium, undergoes multiple scatterings 
and eventually localizes with an exponentially decaying 
density profile [3-5]. The AL is a single-particle inter¬ 
ference effect which cannot be observed directly in solid- 
state systems due to the presence of electron-electron and 
electron-phonon interactions. AL has been widely stud¬ 
ied for various systems including tight binding models 
with diagonal disorder and nearest-neighbor (nn) and/or 
next-nearest-neighbor (nnn) hopping [6-9]. 

Ultracold atomic gases have become a playground 
where complex systems can be simulated and investi¬ 
gated [10-12]. In particular, optical lattice engineering 
allows a high degree of controllability. Various tech¬ 
niques, i.e. periodic modulation of lattice positions and 
amplitudes [13-15] or Feshbach resonance [16], are used 
to effectively tune parameters of a given model [17-21]. 
The level of experimental control and detection allows 
one to build quantum simulators, i.e., experimentally 
controlled systems that are able to mimic other systems 
difficult to investigate directly [22, 23] . Ultracold atomic 
gases are ideal systems for theoretical and experimen¬ 
tal investigation of the Anderson localization of matter 
waves. An experimental observation of AL was first re¬ 
alized in one dimension seven years ago [24] (a closely 
related Aubry-Andre [25] localization was also realized 
[26] independently. For reviews see [27, 28]. Recently AL 
was observed also in three dimensions [29, 30] in speckle 
potentials. In all these experiments, in order to get rid 
of particle interactions, either Feshbach resonances [16] 
were employed or a low atomic density limit was reached. 

Typically in cold atom disorder experiments, as well as 
theoretical propositions, the disorder appears in a diag¬ 
onal form, either on the chemical potential level (lattice 
site energies) like for quasi-periodic lattices [26] or for in¬ 
teractions [31]. Similarly, a diagonal disorder appears for 
a binary disorder that is introduced by interactions with 


other species [32-36]. On the other hand, a more de¬ 
tailed analysis shows that for cold atoms the off-diagonal 
disorder, i.e. a disorder in tunnelings, appears in quite 
a natural way both for the incommensurate superlattice 
potential [37] and for the speckle random potential per¬ 
turbing the optical lattice [38]. However in both these 
cases, the disorder in tunneling is strongly correlated with 
the diagonal disorder. The aim of this paper is to show 
that periodic modulations of interactions allow a trans¬ 
fer of the disorder to the kinetic energy (tunneling terms) 
creating a tunable off-diagonal disorder opening up the 
possibility of its study in controllable cold atoms settings. 

The interest in off-diagonal disorder stems from the 
fact that it can profoundly affect the properties of the 
system. Consider, for example, the case when the disor¬ 
der is purely off-diagonal. There has been a long debate 
about the nature of states in the center of the band. The 
first works [39-41] showed that the localization length 
diverges in the center of the energy band (i.e. E = 0), 
therefore it was argued that a transparency window ap¬ 
pears and even the one-dimensional system exhibits the 
mobility edge. However, in the early eighties it has been 
established that even for a purely off-diagonal disorder 
all states are localized [42-45]. The puzzle of the E = 0 
state was solved showing that the wavefunction envelope 
scales as exp{—j'/N) with N being the system size [42]. 
In the presence of both diagonal and off-diagonal disorder 
all states remain localized [46] except very special cases 
of correlated diagonal and off-diagonal disorder [47] . Ad¬ 
ditional details, especially in the context of conductance 
anomalies, may be found in [48]. 

Possible correlations in the system and/or in the dis¬ 
order may also profoundly affect localization properties 
[49-51]. A famous example is the dimer model [50] in 
which sites may have energies Ca and e;, with the constrain 
that Cfe-sites come only in pairs. This short-range corre¬ 
lation leads to delocalized states. A similar situation oc¬ 
curs for the dual random dimer model (DRDM) in which 
consecutive sites may not have et energy. The model has 
several applications in different areas from DNA studies 
[52, 53] to photonic systems [54] . The cold atom version 
of DRDM has been proposed by Schaff, Akdeniz, and Vi- 
gnolo [55] who have shown that by tuning the interaction 
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between atomic species one may observe the localization- 
delocalization transition and study the extended resonant 
mode. In our proposition we modify the approach pro¬ 
posed in [55] by periodically modulating the interactions 
which allows us a great freedom in changing the relative 
importance of diagonal and off-diagonal disorder in the 
system. For such a model we calculate the approximate 
Anderson localization length and compare it with numer¬ 
ical results. We show that the delocalization window may 
serve as a narrow energy filter for the particles. In addi¬ 
tion we point out that the extended mode is vulnerable to 
effects due to next-nearest neighbor tunneling limiting its 
existence to relatively deep lattices. Details of numerical 
methods needed to calculate the Anderson localization 
length with next-nearest neighbor random tunneling are 
presented in the Appendix B. 


II. THE MODEL 

We consider a simple, standard one-dimensional non¬ 
interacting tight-binding Hamiltonian (with h = \) 

H (etUi - J(ajaj+i -I- h.c.)^ , (1) 

i 

where (a|) denotes an annihilation (creation) opera¬ 
tor of bosons at site i, = a|ai is a particle counting 
operator, are on-site energies while J is the tunneling 
amplitude which will serve as our energy scale. Let us 
apply a periodic modulation of the on-site energies 

Ci —>■ ei(l-I-5sin(a;t)), (2) 

where 6 is the modulation amplitude and w its frequency. 

Hamiltonian (I) has periodic time dependence iL(<) = 
H{t + In such a case we can use the well es¬ 

tablished formalism of Floquet theory [56] (see also [57]) 
extensively used, in the last century, to study the in¬ 
fluence of optical or microwave fields on atoms. Re¬ 
cently it has been applied with great success for con¬ 
trolling properties of ultra-cold atomic systems [13-15]. 
Solutions of equation idt\4’n{t)) = -ff(i)|'0n(i)) have a 
form of Floquet states: \'ipn{t)) = e“*®”‘|u„(t)), where 
is called the quasienergy and |it„(t)) have period¬ 
icity of the Hamiltonian. The Floquet theorem is a 
time analogue of the Bloch theorem for spatially peri¬ 
odic potentials. Although we can not treat \tj}(t)) as 
eigenstates of H{t), \u{t)) are eigenstates of the Flo¬ 
quet Hamiltonian 'H{t) = H{t) — idt existing in the ex¬ 
tended space of T-periodic functions. In that space we 
can number states using a new quantum number m £ Z: 
Wri) = where |u™) is the eigenstate to the 

eigenenergy e™ = e„ -bwm. This whole class corresponds 
to one physical state \ipnit)) as adding w to in phys¬ 
ical space only takes us to the next “Brillouin zone” for 
quasienergies. Therefore, it is sufficient to find a block- 
diagonal form of a Hamiltonian (in m-ordered basis) and 


consider only one block corresponding to a single “Bril¬ 
louin zone”. Unfortunately, couplings between different 
blocks make the block-diagonalization a formidable task. 
One may, however, make a unitary transformation U: 


n' = unu\ 


U = exp 


,cos(a;t) 
—id - 

LJ 



( 3 ) 


T-L'{t) = ^ [eiTii - (e 




5 cos(u)t) 


ajai+i -I- h.c. 


( 4 ) 

In effect, for frequencies w ^ 1 (in units of J) one we can 
neglect the couplings between different Floquet blocks 
and consider only one diagonal block governing the long 
term (t ^ 1/w) dynamics: 


HeS — rji 


dt'H'it) = ^ - (Lafoi+i -|-h.c.)^ 


( 5 ) 


where 

ti = Jo (ci+i — Ci)^ (6) 

is the effective position dependent hopping and Jq is the 
zero-th order Bessel function. It has been verified experi¬ 
mentally [14, 15] that the effective time-averaged Hamil¬ 
tonian governs the dynamics of periodically driven sys¬ 
tems for long times. In particular, localization properties 
of eigenstates of (5) will be shared by the Floquet eigen¬ 
states of the original Hamiltonian (1) [with the fast mod¬ 
ulation (2)[ after averaging them over the period. This 
may be understood by inspecting the unitary operator U 
in (3) and observing that it consists of a product of local 
operators acting on single sites. Thus U adds only local 
phases to single particle states - that cannot affect the 
resulting density distribution. Further detailed discus¬ 
sion of the construction of the effective Hamiltonians for 
high frequency periodically driven systems may be found 
in [58-60]. 

Note that for a uniform system with all being the 
same, the tunneling, (6) is unaffected, while ti changes if 
on-site energies vary from site to site. Such is the case 
for super lattices [61, 62] or external potentials such as a 
harmonic trap or a linear tilt [13]. In this work we will 
consider the on-site energy variations due to disorder. 

It is worth noting that it is usually possible to adi- 
abatically pass from eigenstates of one ideff to another 
(e.g. for different 6) if the change of parameters is slow 
enough [63]. Thus it is possible to prepare a time in¬ 
dependent system and turn on periodic modulation, or 
change modulation parameters during experiments. For 
completeness let us note that, for the tight-binding de¬ 
scription in terms of a single lowest band to be valid in 
(1-5), w, while much larger than the tunneling amplitude, 
must be smaller than the energy separation to the excited 
band [13, 64] . 
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To create disorder we consider two species of atoms re¬ 
pulsively interacting with each other but allow only one 
of them to move freely on the lattice. The frozen atoms 
(denoted with a superscript /) give rise to a binary disor¬ 
der [33-35, 65]. The mobile particles are assumed not to 
interact among themselves. They can be spin-polarized 
fermions or bosons with interactions tuned off by mi¬ 
crowave or optical Feshbach resonance [16]. The dynam¬ 
ics of the system can be described by a single particle 
Hamiltonian with the on-site energy Ci = Vnj, where V 
denotes an interspecies contact interaction. When frozen 
particles are fermions or strongly repelling (hard-core) 
bosons their occupation per lattice site nj takes either a 
value of zero or one and the on-site energy takes only two 
values Ei G {0, F}. Further, we consider the particular 
case of DRDM when two adjacent sites cannot be occu¬ 
pied by frozen atoms simultaneously [50] . For cold atoms 
DRDM can be created by the method described in [55]. 
Referring the reader to that paper for details, let us men¬ 
tion only, that to assure no heavy particles reside in the 
neighboring sites one may first trap the heavy species in 
an auxiliary long-wavelength lattice with the lattice con¬ 
stant being e.g. three times bigger than the final lattice 
to be considered. If these heavy particles are strongly re¬ 
pelling and for sufficiently low densities one may assure 
no double occupancy. Only then one can switch on the 
desired shorter wavelength lattice which holds the mobile 
component. 

In [ 66 ] the authors considered fast periodic modulation 
of interactions induced by an appropriate periodic mod¬ 
ulation of the magnetic field B{t) = B{t + T) with period 
T = ^ for interacting bosons in the optical lattice. We 
assume a similar mechanism applied to light-heavy par¬ 
ticle interactions resulting in periodic (for simplicity as¬ 
sumed to be harmonic) time dependence of site energies 
in the form ei{t) = n{(V q + Fisin(a;t)), where Vi is an 
amplitude modulation strength and uj is the modulation 
frequency. For magnetic field values close to the Fesh¬ 
bach resonance [16] large variations of interactions may 
be expected. In particular choosing the value of the mean 
magnetic field around which the oscillations occur one 
may vary at will the relative importance of the Vq and 
Vi coefficients. Note that we can adjust the magnetic 
field for this purpose since we assume that the mobile 
particles are either fermions or bosons with interactions 
turned off by either microwave or optical Feshbach reso¬ 
nance [16]. 

Compared to the general model discussed above we 
have Ei = n{Vo and S = Fi/Vb- As n{ G {0,1}, the 
effective tunneling, after time averaging, can only take 
two values: 

^fO, ifnf =0 1 , ifnf = nf+i 

\Fo, ifnf = l’ ' ll- 7 , 

(7) 

where 7 = 1 — Jo (hi/w) measures the strength of 
off-diagonal disorder and varies in range from zero to 
about 1.4 (to the minimum of the Bessel function, 1 -|- 


Tah\^Jo{x)). 


III. ANDERSON LOCALIZATION LENGTH 


To calculate Anderson localization length let us start 
with the time-independent Schrodinger equation for the 
disordered tight binding Hamiltonian (5) 

- U-itfi-i + Eilfi = Eifi, ( 8 ) 


where E is the eigenenergy. In the regime of small di¬ 
agonal and off-diagonal disorder (Vq 1,7 <C 1) we 
can assume that E is approximately given by the disper¬ 
sion relation E Ri —2 cos(fc), where k is the quasimomen¬ 
tum. In order to define Anderson localization length in 
a system with ofl-diagonal disorder we utilize a unitary 
transformation tpi = (firji, where = l/[ti? 7 i_i]. We 
transform equation ( 8 ) to diagonal form: 

+ f’i-l + = Oj (9) 


where Vi = |? 7 iP(ei -I- 2cos(A:)) is a new effective diagonal 
disorder. In the considered DRDM model rji = 1/ti. 
Equation (9) can be expressed as a two dimensional 
Hamiltonian map with position and momentum of the 
form Xi = 4>i, Pi = {(fi cos(fc) -|- (fi-i)/ sin(fc): 


(xi+i = -{pi + AiXi)sm{k) + XiCOs{k) 

= {pi + AiXi)cos{k) + XiSm{k), 

where 


Ej - 2(1 - |tip)cos(fc) 
sin(fc) 


( 10 ) 


( 11 ) 


The map (10) describes the behavior of a harmonic oscil¬ 
lator under periodic delta kicks with amplitude Ai. Ex¬ 
pressing the map in action-angle variables and iterating 
it one may estimate the localization length in the limit 
of small diagonal and off-diagonal disorder. The details 
are given in the Appendix. We obtain the approximate 
inverse localization length 

^-l P (Fq- f 27(2 - 7) cos(fc))^ 

H + pY 8sin2(fc) 

/ ^ p(p + cos( 2 fc)) \ 

V l + p 2 + 2 pcos( 2 fc)y ’ 

where p = p/{l — p) and p is the mean occupation number 
of frozen particles. A simple analysis of equation (12) in¬ 
dicates an anomalous behavior of Anderson localization 
length due to the existence of a transparent mode with 
momentum kf given by condition cos(fct) = — 
which the localization length diverges. Frozen particles 
in a lattice increase the on-site energy and a motion of 
mobile particles can be seen as a sequence of scattering 
processes on potential barriers. In the DRDM two frozen 
particles are separated by at least one lattice site and the 






resonance is caused by a lossless transmission through a 
single barrier. The existence of the transparent mode 
is given by the condition Vb < 27(2 — 7 ) and the wave 
with momentum kt travels through the sample without 
reflection. A similar anomalous mode was observed in 
Ref. [55]. The on-site energy modulation allows us to 
change the off-diagonal disorder significantly while keep¬ 
ing the amplitude of the diagonal disorder small (con¬ 
trary to models deriving the changes of the tunneling 
from changes of an effective lattice shape only [55, 67]). 
In Fig. 1 we present the numerically calculated Anderson 
localization lengths (using the standard transfer matrix 
method, see e.g. [ 68 ]) and compare them with the ana¬ 
lytical expression (12). The left panel corresponds to the 
weak diagonal (Vb € {0.05,0.1}) and weak diagonal and 
off-diagonal (7 € {0.05,0.1}) regime. The occupation 
of frozen particles is equal to p = 1/3. The theoreti¬ 
cal predictions are in good agreement with the numerical 
calculations. The right panel corresponds to the weak 
diagonal (Vb = 0 . 1 ) and strong off-diagonal (7 = 0.9) 
regime. The width of the divergent Anderson localiza¬ 
tion length window is significantly narrower than for the 
weak disorder case. The position of the transparent mode 
kt is properly described by Eq. (12) while the Anderson 
localization length shape is not. Due to the existence of 
the transparent mode in the system, one can expect that 
the evolution of an initial wave packet with a given mo¬ 
mentum distribution results in Anderson localization of 
all momenta except those very close to the transparent 
mode kt- Disorder effectively works as a band-pass fil¬ 
ter for momenta. To verify this behavior we integrated 
the Schrodinger equation in time for the Hamiltonian (5) 
with the initial state consisting of one particle in the cen¬ 
ter of the lattice with uniform momentum distribution. 
Fig. 2 presents momentum distributions of Anderson lo¬ 
calized atoms and those that escaped from the disorder 
after time evolution t = 1000. The momenta of escaped 
atoms reveal narrow peaks at the position of the trans¬ 
parent mode while the momentum distribution of Ander¬ 
son localized atoms reveal dips, positions of which agree 
with the former peaks. The lattice has N = 1024 sites 
with filling p = 1/3. We choose two sets of parameters 
for the amplitude of diagonal and off-diagonal disorder 
(Vb = 0 . 1,7 = 0.5) and (Vb = 0 . 5,7 = 0.3). For such 
parameters the transparent mode appears for kt ~ 1.64 
and kt ~ 2.08 respectively. Our system presents, there¬ 
fore, a promising candidate for obtaining a highly con¬ 
trollable monochromatic gun for mater waves. A similar 
mechanism was observed in [69] where coherent multiple 
scattering processes determined the emitted matter-wave 
mode. In the following we shall inspect residual effects 
that may affect the existence of the transparent mode for 
a realistic system. In particular, in the next Section we 
analyze the influence of the next-nearest neighbor tun¬ 
neling on Anderson localization length in the regime of 
strong off-diagonal disorder. 
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FIG. 1. (color online) Anderson localization length obtained 
from transfer matrix calculations (solid lines) and from ( 12 ) 
(dashed lines). Left panel: Anderson localization length as a 
function of quasimomentum for two sets of diagonal and off- 
diagonal disorder amplitude: {Vb = 0.05, 7 = 0.05} (black) 
and {Vb = 0 . 1,7 = 0.1} (red). Filling p = 1/3. Right panel: 
Anderson localization length as a function of quasimomentum 
from transfer matrix calculation (black solid line) for diagonal 
disorder amplitude Vb = 0 . 1 , off-diagonal disorder amplitude 
7 = 0.9 and occupation p = 1/3. For comparison with (12) 
we present analytical expression (black dashed line). 



quasimomentum 

FIG. 2. (Golor online) Horizontal dotted thick line repre¬ 
sents the initial wave packet momentum distribution. Solid 
lines present momentum distribution of wave function out¬ 
side area of the disorder while dashed line the distribution of 
momenta that remain in the system. N = 1000 lattice sites 
are considered with filling p = 1/3. Two cases are plotted. 
The data for diagonal and off-diagonal amplitudes Vb = 0.1, 
7 = 0.5 are represented by red (thick) curves while those for 
Vb = 0.5, 7 = 0.3 by black (thin) lines. The evolution time 
is t = 1000. The position of dips in momentum space agree 
with position of transparent mode kt « 1.64 and kt « 2.08 
respectively. One can observe two peaks in the momentum 
distributions of escaped atoms. Due to a marginal difference 
in the dispersion relation between in and outside of the dis¬ 
order, the positions of peaks and dips are slightly shifted. 
These results are obtained as the average over 1000 disorder 
realizations. 
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FIG. 3. Panel (a): One dimensional gas on a lattice with 
nearest neighbor hopping (solid line) and next-nearest neigh¬ 
bor hopping (dotted line). Panel (b): The same system can 
be viewed and described differently, namely as a stripe with 
non uniform tunneling. The advantage is that the hopping of 
particles takes place between neighboring slices only. 
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FIG. 4. (color online) Anderson localization length as a 
function of quasimomentum for diagonal disorder Vo = 1.95 
(top panel), Vb = 1 (right bottom) and Vo =0.1 (left bot¬ 
tom), 7 = 0.9, occupation p — 1/3. The value of the 
next nearest neighbor hopping is t' = 0.01 (black solid line), 
t' = 0.005 (blue dashed line), t' = 0.0025 (red dotted line) and 
t' = 0.0005 (black dashed-dotted line). The resonance dimin¬ 
ishes when t' grows. The effect is the strongest for Vo = 1.95 
because we are close to a regime without a delocalized mode. 


the Hamiltonian (1) 

H —> H 1' ^^(a|ai -)-2 + h.c.). (13) 

i 

Under the same unitary transformation (3) and after time 
averaging we obtain the effective Hamiltonian in the form 
of (5) plus the additional term reading 

Y^t'i{alai +2 + h.c.) (14) 

i 

with t' = t'Jo i^i +2 — £i)) ■ We numerically calcu¬ 
late Anderson localization length for this new effective 
Hamiltonian. The method we use is applicable to finite 
but arbitrarily long-range hopping. The transfer matrix 
approach is not directly applicable so we transform the 
problem from a one dimensional to a two dimensional 
stripe. This could be achieved with a simple winding of 
a string, see Fig. 3. In a stripe all particle hoppings are 
between the neighboring slices, which simplifies numeri¬ 
cal calculations. 

In the Appendix we present the derivation of the equa¬ 
tion for a Green’s function between the first and the last 
slice of the stripe. Knowing the Green’s function we can 
easily calculate the localization length [70, 71]. 

As exemplary parameters for numerical analysis we 
choose Vo = {1.95,1,0.1} and 7 = 0.9 in order to obtain 
a regime where a next-nearest-neighbor hopping signif¬ 
icantly affects the dynamics of the system. The results 
are presented in Fig. 4. We observe that even extremely 
small values of t' significantly affect the resonance. The 
presence of a small but non-zero next-nearest-neighbor 
tunneling t' significantly reduces the localization length 
for the transparent mode kt due to a non-zero proba¬ 
bility that a mobile particle scatters back on a frozen 
particle through a next-nearest-neighbor hopping. This 
process can appear when two frozen particles are sep¬ 
arated by a single site. Indeed, when we exclude such 
configurations, the resonance reappears. In deep lattices 
next-nearest-neighbor tunneling also reduces the local¬ 
ization length for the transparent mode, but the effect is 
weaker and may not affect the selective emission due to 
a finite optical lattice size (when this size is smaller than 
the localization length). 


V. CONCLUSIONS 


IV. NEXT-NEAREST NEIGHBOR 
TUNNELINGS 

For lattice depths U ~ 10-20 K/j, where Er is the 
lattice recoil energy and the typical optical laser wave¬ 
lengths, next-nearest neighbor hopping amplitude t' 
ranges between 0 . 01 - 0 . 001 , as can easily be estimated 
with appropriate Wannier functions. This is why in typi¬ 
cal situations long range hopping is often negligible. Still 
let us add the next-nearest neighbor tunneling term to 


In this paper we proposed a method to realize a con¬ 
trollable off-diagonal disorder with binary random poten¬ 
tial resulting from time-periodic modulation of interac¬ 
tions between mobile and frozen particles in an optical 
lattice. Since no interactions between mobile particles 
are taken into account they may be assumed to be spin¬ 
less (spin-polarized) fermions. One could also imagine a 
scheme with bosons with interactions tuned off by some 
microwave or optical Feshbach resonance [16] (note that 
we already assume using a magnetic field for a standard 
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Feshbach tuning of different species interactions so this 
method cannot be used simultaneously to control light- 
light particle collisions). The presented method allows us 
to obtain models with strong off-diagonal and weak diag¬ 
onal disorder in a broad regime of relative (off)-diagonal 
disorder amplitudes. In a regime of weak off-diagonal 
disorder, an analytical expression for the Anderson lo¬ 
calization length is in very good agreement with numeri¬ 
cal simulations. Moreover, we indicate how DRDM with 
random hopping can work as a tunable band-pass filter 
for matter-waves. The momentum of escaped atoms re¬ 
veal narrow peaks in the position of the transparent mode 
in momentum space while momentum distribution of the 
Anderson localized atoms reveals dips whose positions 
agree with former peaks. We indicate the importance 
of the next-nearest-neighbor hopping for the localization 


length of the resonant extended mode appearing within 
dual dimer random model for strong diagonal disorder. 
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VII. APPENDIX A:HAMILTONIAN 
APPROACH TO ANDERSON LOCALIZATION 
LENGTH 


To find the localization length starting from the map 
(10) [47, 72, 73] it is more convenient to express it in 
the action-angle variables (r, 9) using the transformation 


X = rsin(0), p = rcos(0): 

s'm{6i+i) = D~^{sm{6i - k) - Ai sm{di) sin(fc)) 
cos( 0 i+i) = D~^{cos{9i — k) -i- Ai sin( 0 i) cos(fc)) ’ 

(15) 


where 


D, = 


D-ri 


= JI ^-2AiSin{29i) + A^siT?{9i). (16) 


We define Anderson localization length as 

N 


A lim — y^lnUi, 
Af^-oo N ^ 


(17) 


i=l 


and after the expanding logarithm to the second order in 
Ai we get 

A-' = ^(Al?) + i(A,sin(20O), (18) 

where (...) stands for averaging over i. In order to cal¬ 
culate the ’kick-angle’ correlator {Ai sin( 20 i)) we expand 
the map (15) to second order in 0^: 


— 9i-\ — fed- Ai- 


i-l) 


sin^(0; 

^ sin(fc) 


(19) 


and express {Ai sin( 20 i)) in terms of the preceding 
(Ai sin(20i_i)) ’kick-angle’ correlator. Let us introduce 
the correlation of the kick’s strength Ai with angle 9i 


2i 

W) 




( 20 ) 


where {Af) is variance of the Ai. Multiple applications 
of eq. (19) to a„ give us the recursive relation 


_ 2jh 

O^n—l — 6 H" 


( 21 ) 


where qn = {AiAi-n)/{Af) is the autocorrelation of Ai. 
From the definition of a„ we can notice that 

{A, sin(20,)) = 5R , (22) 

where from ( 21 ) we obtain ag = 

The inverse Anderson localization length takes the 
form [72] 


A 


-1 


= + 2 ^ qn(k) cos( 2 kn) 


(23) 


In the specific case of DRDM , the variance of Ai and 
the autocorrelation of Ai read, respectively 


{Ai) = 


p {Vq -b 27(2 - 7) cos(A:))^ 


(1 + p)^ 

= (-!)>" 


sim(fc) 


(24) 


and the approximate expression for the inverse local¬ 
ization length takes the form ( 12 ). 









VIII. APPENDIX B:ANDERSON 
LOCALIZATION LENGTH FOR HAMILTONIAN 
WITH LONG-RANGE RANDOM HOPPING 

In this section we describe a method of calculating lo¬ 
calization length in a one dimensional system with long- 
range tunneling. The idea is to transform the problem 
from a one dimensional to a two dimensional stripe. It 
could be achieved with a simple winding of a string, 
see Fig. 3. 

A. Mapping to a two-dimensional stripe 

To start let us consider a one dimensional gas on a 
string of length N. We allow long-range tunneling up to 
the M-th neighboring site. This geometry is equivalent 
to a two dimensional stripe of length N = N/M and 
width M. The evolution of particles in such a stripe is 
characterized by a Hamiltonian 

N N-1 

= ^ (h„ + /i.c.), (25) 

n—1 n—1 

where 

M 

^ e^^'>\n,Tn){n,Tn\ 

m—1 

M-1 

~ +^-C-) ( 26 ) 

m—1 

is a standard tight binding Hamiltonian for the n-th slice 
and 

M M' 

Vn =+ (27) 

m—1 m' — l 

describes particles hopping onto the n-th slice. 
denotes tunneling between |n,m) and |n',rn'). Suppose 
now that we add one extra slice to the system. A new 
Hamiltonian is straightforward 

HiN + l) = HiN)+HN+i + VN + Vl. (28) 
For later convenience, let us define operators 

H°{N + 1) = H{N) + Hn+i (29) 

ViN) = Vn + V^, (30) 

and finally the Hamiltonian reads 

H{N + 1) = H°{N + 1) + ViN). (31) 

B. A recursive equation for a Green’s function 

Our goal is to find the Green’s function + I) for 

the Hamiltonian H{N+1). By definition, a Green’s func¬ 
tion can be obtained from a resolvent of a Hamiltonian 


G^(IV-|-1) — ^Ih^ G£;+ie(V-|-1), (32) 

where a resolvent satisfies 

(z-iJ(V + I))G,(iV + l) = 1. (33) 

Let’s start with an equation for a resolvent of H{N+1): 

G,{N+1) = G°{N+1)+G°{N+1)V{N)G,{N+1), (34) 

with G° being a resolvent of H^. The resolvent equa¬ 
tion (34) can be derived from a simple operator identity 

with a substitution 

A = z- H{N + 1), B = z- H°iN + 1). (36) 

Since we are interested in transport properties of the 
system, we need to know the Green’s function matrix ele¬ 
ments between the first and the last slice only. Therefore, 
we want to obtain (I|G 2 (iV -|- 1)|V-|- I). To simplify the 
equations let’s introduce a notation: 

Gn,m = (n|G^(m)|m). (37) 

From equation (34) we get 

Gi,n+i = Gi,nVn Gn+i,n+i, (38) 

where we used an observation that 

G°(iV + l) = G,(iV) +- 1 -, (39) 

z - Hn+1 

which stems from the fact that H{N) and Hn+i act in 
orthogonal subspaces of the Hilbert space, see (29). 

In order to get a recursive equation for Gi^at, we need 
to calculate Gat+i^at+i matrix. It can be obtained from 
equation (33) by multiplying it from the right side by a 
projection P = \N + I)(V-|- 1|, and by P (or Q = 1 — P) 
from the left side: 

P{z - HiN + 1))(P + Q)G,(iV + I)P = P, 
Q{z-H{N + l)){P + Q)G4N + l)P = 0. (40) 

Solving a set of equations (40) one gets 

GAr-|-i,Ar-|-l = -;-• (41) 

Z — Pat- 1-1 — bA)-GAr,AfVAr 

Gombining equations (38) and (41) 

Gi,Ar-i-i (^z — Hn + 1 — V^Gn,nVn'^ = Gi^nVn (42) 

and, again from (38), extracting Gn,n we finally obtain 
a recursive equation: 

An+2 = {E- Pjv+OV^^Aat+i - (43) 
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with An = ]si-2 ■ initial values of Ajq are not 

relevant, but it is convenient to choose them as 


^0 — 0, Ax — Vq. 


(44) 


In equation (43) there is no singularity, hence we could 
replace z with E. From this equation we can calculate 
Gl.N^ which is connected with the localization length of 
the system. 


C. Calculation of the Anderson localization length 

The Anderson localization length Am in a stripe of 
width M is defined as 


- lim ilnTr|C?i„p. 
Am n^oo n 


(45) 


Since our system is in fact one dimensional the local¬ 
ization length should not depend on the length of the 
stripe. When a particle reached the n-th slice it covered 
n = M ■ n lattice sites. Hence, the localization length A 
in one dimension equals A = MXm- 

The localization length can be obtained from the re¬ 
cursive equation (43). We solve the equation iteratively. 
The problem with this equation is that the elements of 
An grow exponentially for large n, so it requires some 
regularization. Therefore, in each step we multiply the 
both sides of the equation by some matrix i?„. Starting 
from n = 1 and defining = AkRi 


A 3 = (A - H2)VG^A2 - HiVp-iAi 


X Ri 


= [E- H2)Vx~^Ai^'> - 


(46) 


To avoid the exponential growth we put = A 2 ^ and 
a('UA3A2-\ aW = 1, aW=AiA 2 -\ (47) 
We repeat this procedure in every iteration so that 




Rm with Rn = 


^n+1 


satisfies 


Ai% = {E- Hn+x)V-^ - 


(48) 


^n+2 

Let US also define a matrix 

bM = bn = (49) 


where 11-11 = V^Tr| • p is a matrix norm. matrix 

turns out to be very useful, because 


bn = \\B^^-^^Rn\\ = - - 

On-1 


bn — 1 




a: 


(n-2) 

n-t-1 


-1 


1 


bn—lbn—2 
Continuing this, we notice that 


^(n-3) 


4("-3) 

2^71-1-1 


-1 


(50) 


\\A-l,\\=bxb2...bn, (51) 

lnTr|Gi_„p = 2 (ln 6 „+i -|- ... - 1 - In 61 ). (52) 

The Anderson localization length can be expressed as: 
n 


X = —M lim 


"j E bn.bn-\-l- (^^) 


n-s-oo c„+l 

Iterations should be continued until A converges within 
a desired precision. 


























